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Abstract 

> 

■ I outline various derivations of the non-Abelian Kubo equation, which gov- 
CN , erns the response of a quark-gluon plasma to hard thermal perturbations. In 

the static case, it is proven that gauge theories do not support hard thermal 
<^ • solitons. Explicit solutions are constructed within an SU(2) Ansatz and they 

are shown to support the general result. The time-dependent problem, i.e., 
non-Abelian plasma waves, has not been completely solved. We express and 
motivate the hope that the intimate relations linking the gauge-invariance 
condition for hard thermal loops to the equation of motion for T = 0, topo- 
logical Chern-Simons theory may yield new insight into this field. 
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1. The Quark- Gluon Plasma 



The study of the quark-gluon plasma (QGP) phase of high-temperature QCD 4 is rele- 
vant both to to big-bang cosmology and to heavy ions collisions. Experiments at CERN and 
Brookhaven have provided compelling evidence for the existence of an intermediate "phase" of 
QCD during such collisions, with an energy density ~ 1 GeV/fm 3 (see for instance It is 
however not clear if this "phase" is a thermodynamically distinct one. 

A heavy-ions interaction, after the primary collision of ions bunches, forms a thermalized 
state which lasts for most of the total interaction time. It is in this thermalized state that the 
QGP "phase" is thought to be formed. The interaction ends by hadronization of the thermalized 
constituents. Each step of the process can be tested specifically [l], However, the methods 

1 Talk given at THERM095, the 4th International Workshop on Thermal Field Theories and Their Appli- 
cations, Dalian (P.R. China), August 6-10, 1995. This work has been supported in part by the Swiss National 
Science Foundation. 
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available for analysis are empirical, and one feels compelled to work towards developing a more 
satisfactory theoretical framework for the investigation of the QGP. 

2. High- Temperature Perturbative QCD 

High-temperature perturbative QCD provides a natural setting for the theoretical study 
of the QGP. We shall adopt the hypotheses that (i) QCD is in its deconfmed phase and that (ii) 
the plasma is collisionless, so that it can be regarded as a dilute gas of quarks and gluons. The 
dominant interaction is between hard colored plasma particles and the soft mean color field. 
"Hard" denotes an energy scale of the order of the temperature, e ~ T, while "soft" refers to 
the energy scale e ~ gT, g being the gauge coupling constant. It is to expect that motion on 
a distance scale ~ 1/gT, the inverse soft energy scale, should involve coherently many hard 
particles, giving rise to color polarization effects in the QGP. Indeed, the real part of the color 
polarization tensor describes Debye screening of the color charges, while its imaginary part 
accounts for Landau damping of the mean color field. 

If zero-temperature perturbative QCD is naively extended to the T > case, physical 
quantities turn out to be gauge-dependent. This old puzzle has been solved by the discovery 
of the so-called "hard thermal loops" (HTLs) phenomenon ||: at high temperature, effects 
of leading order in g can arise at any order of perturbation theory, i.e., loop diagrams exist 
that are comparable to the corresponding tree amplitude. Such diagrams contain one-loop 
integrals with soft external lines and hard internal (loop) momenta, and are termed "hard 
thermal loops". In order to account consistently for the orders in g, one is led to resum the 
perturbative expansion. 

Hard thermal loops represent thermal, not quantum, corrections. They describe polariza- 
tion effects in the plasma, as well as the propagation of plasma waves, plasma solitons, etc. The 
hard thermal loops are the dominant diagrams at high temperature. They are proportional to 
the Debye mass squared m D = {N + ^f~), where N is the number of colors and Np the 
number of flavors. Note that the Debye mass is proportional to the soft scale gT; mo represents 
the smallest possible frequency for time- dependent (gluon) propagation. The Debye screening 
length is just its inverse, Xd = ™~d ~ ^/gT. 

3. The Non-Abelian Kubo Equation 

Hard thermal loops are generated by an effective action ||, |]] : 



The first term in the right side describes color-electric screening, and W / (A + ) is a non-local 
functional of A + = ^=(l;q) • A, where q is an arbitrary unit 3-vector, q 2 = 1. The measure 
denotes integration over the solid angle spanned by q. 

The response of the QGP to an external perturbation is governed by a Yang-Mills field 




(3.1) 
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equation 
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Dfj,F^ — ^ induced > (3-2) 
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with the induced hard thermal current J^ duced = — ^-Fhtl(^4) as its source^. In order to 
close the system of differential equations for A^, we need a constraint on the current. Exploiting 
the analogy between the gauge invariance of HTLs and Chern-Simons theory^ || , one derives 
the induced current in the form || 



^induced — J ^ \y 



a_(x) — A-{x) 



Qt 



a + (x) — A + (x) 



(3.3) 



where Q± = -^=(1, ±q), q 2 = 1, are light-like 4- vectors, A± = Q±A^ are light-like projections 
of Ap, and a± are constrained, introducing d± = Q±d^, by 

d + a_ - d_A + + [Ah o_] = , 9+A_ - <9_a+ - [A_, o + ] = . (3.4) 



From now on, we shall call non-Abelian Kubo equation the system (|3.2|) - (|3.4j) , which governs 
hard thermal processes in the QGP. 

Alternative derivations of the non-Abelian Kubo equations have been given. 

1. The non-Abelian Kubo equation can be derived from quantum kinetic theory 0. In 
that approach, one starts from the Schwinger-Dyson equations, and uses a kinematic approxi- 
mation scheme in which derivatives with respect to center-of-mass and to relative coordinates 
carry different (/-dependences. From the resulting quantum kinetic equations, one is able to 
derive F H tl(^4), and hence an expression for the induced current. 

2. The non-Abelian Kubo equation can also be derived, as done in ||, by applying the 
action principle to the Cornwall- Jackiw- To mboulis composite effective action T[<f), G^] |§, a 
functional of a generic field <fi (which stands for gluons, quarks, antiquarks, ghosts, etc) and 
of the two-point functions G^. Requiring that r[</>, G^] be stationary with respect to <f>(x) 



and G<j,(x,y) yields, respectively, the equation ( |3.2j ) and a set of conditions on the two-point 



functions which can be cast in the form of a constraint on the induced current 

This approach is close to the one of 0, of which it adopts numerous steps, including its 
kinematic approximation scheme. The composite action approach presents the advantage of 
providing an action principle for hard thermal processes in the quark-gluon plasma. Whether, 
and which, physical hypotheses can be expressed at the level of the composite action functional 
is an open question. 

3. More recently, the non-Abelian Kubo equation has been derived from classical trans- 



port theory [10|. This is the subject of the talk by C. Manuel at this meeting. One starts from 



classical particle dynamics as described by the Wong equations; the color charge is considered as 



1 Note that (|3.2|) is just the equation of motion for the gauge field, when one adds the zero-temperature QCD 



action to the HTLs generating functional, i.e., (3.2) is equivalent to requiring xa~(^QCD + Thtl) 



0. 



See Section 5 below for a more detailed account of this analogy. 
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a classical phase-space variable. Transport theory yields the non-Abelian Vlasov equations, and 
the induced current is expressed in terms of the one-particle distribution functions. Expanding 
these distribution functions in powers of the coupling constant g yields, at leading order, the 
gauge invariance condition for the generating functional of hard thermal loops Thtl(^4) ( |3.1|) . 
Subsequent developments that lead to the form (|3.3| ) -(|3.4| ) for the induced current, and hence 
to the non-Abelian Kubo equation, are described in ||. 

4. Solutions of the Non-Abelian Kubo Equation: QGP 
Response 

It is of obvious interest to discuss solutions of the non-Abelian Kubo equation (|3.2|)-( |3~4]) . 
In the Abelian, electrodynamical case this is easy to do, since the conditions fl3.4|) can be readily 
solved for a±, and the solutions of the linear problem are the well-known Abelian plasma waves 



TT| . The non-linear problem of finding non-Abelian plasma waves is much more formidable. 



Also, one inquires whether the non-linear equations support soliton (and instanton) solutions. 

Let us first consider the static case. Acting on time- independent fields A± = A±(q), the 
derivatives in ( |3.4| ) become d± = ±-^=q_- V = ±<9 q ; therefore the constraints Q3.4Q can be written 

as 

d (l {A + + a_) + [A + ,A + + a„) = , <9 q (A_ + a+) - [A_, A„ + a+) = . (4.1) 
These constraints are solved by a_ = — A + , resp. a + = -A-. (U) then yields 

^duced(q) = - / ^ m + Q-) m + Q V MM (4.2) 

for the static inducent current. With Q + + Q = and Q° + + Q°_ = y/2, one computes 
•induced = ( — 2 A ; 0). The response equations ( |3.2| ) become, in the static limit: 

DiE? = -m 2 D A , e l3k D 3 B k - [A , £*]=(), (4.3) 

where E l = F l ° and F u = —e^ B , and the Debye mass mo plays the role of a gauge invariant, 
color-electric screening mass. 

In the study of static excitations, one is of course interested in the issue of possible (static) 
soliton solutions. It has been established in || that the non-Abelian Kubo equation does not 
possess any finite energy static solutions. The proof goes as follows. Consider the symmetric 
tensor 

9 ij = 2 Tr ( y E i E j + B i B j - y (E 2 + B 2 + m 2 D A 2 )^j . (4.4) 
Using (|4.3|) one verifies that for static fields dj O^ 1 = 0. Therefore 

(PrO™ = f <i''r<)j{r' ()'' ) = f dS'xW . (4.5) 



Moreover, the energy of a massive gauge field (with no mass for the spatial components) can 
be written as 



S = I d*r\ -Tr [E 2 + B 2 + ^(A^)^J + Tr (m D A + -±- ) \. ( I.G) 
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The second trace in the integrand enforces the first of constraints (|4.3| ). Consequently, on the 
constrained surface the energy is a sum of positive terms, 

8 = Jd 3 r {-Tr (e 2 + B 2 + m 2 D A 2 )} , (4.7) 

and E, B and A must decrease at large distances sufficiently rapidly so that each of them is 
square integrable. This in turn ensures that the surface integral at infinity in (|4.5| ) vanishes, 
so that static solutions require / d 3 r 8 U = 0. On the other hand, from (O), we see that 6 U is 



a sum of positive terms, 8 n = — Tr (E 2 + B 2 + ?>m 2 D A^), hence we must have E = B = and 
Aq = 0, i.e., there are no (non-trivial) finite-energy solutions to ( [4 . 3| ) . A similar argument H] 
shows the absence of "static" instantons. 

I want to analyze now in some details possible solutions of the static response equations, 
and present numerical results || that support the general proof of absence of (static) solitons. I 
consider for simplicity the radially symmetric restriction of the static response equations 
in the SU (2) case. Radially symmetric SU (2) gauge potentials have the formP] 

A a = f a girl , Af — r - rv) + e« f j l ~^ r) . (4. 



We substitute this Ansatz into ( |4.3|) . The resulting equations give us the freedom to set (f) 2 to 
zero. Rescaling x = mr and defining J(x) = g(r), K(x) = <pi(r), we obtain a system of coupled 
second-order differential equations 

x *f- J = (x 2 + IK 2 ) J , x 2 ^- 2 K = (K 2 - J 2 - 1) K . (4.9) 

dx z ax z 

These equations possess two exact solutions: 

(1) the Yang-Mills vacuum, with J = 0, K = ±1, and 

(2) the Wu-Yang monopole, plus a screened electric field: J = Joe~ x , K = 0. 

Let us investigate the asymptotic behaviour of the system (|4.9p . At x — > oo, the regular 
solution tends to the Yang-Mills vacuum, with J approaching its asymptote exponentially^. 
Near the origin, i.e., at x — > 0, J and K behave either like the Yang-Mills vacuum, or approach 
the Wu-Yang monopole, as 



J(x)->J + ..., K(x) ^ K ^cos^Ji + ^ In^j + ... . (4.10) 

Only the vacuum alternative at the origin could lead to finite energy. However, since 
we must choose one of two possible solutions at infinity (obviously we pick the regular one), 
the behavior at the origin is determined and can be exhibited explicitly by integrating the 
equations ( |4.9| ) numerically. Starting with regular boundary conditions at infinity, we find that 
the monopole solution is reached at the origin, with K vanishing as in (|4.10p ; see Figure 1 

3 A residual gauge freedom has been used to eliminate a term proportional to f a f l . 
4 There is also a solution with J growing exponentially, which we do not consider. 
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below. This result is consistent with the analytic proof |§ that there are no finite energy static 
solutions in hard thermal gauge theories. 

Much less is known about time-dependent solutions of the non-Abelian Kubo equation 
(|3.2|) - ([57^ ). The special case of time-dependent, space-independent solutions has been investi- 
gated in The Ansatz A^(x) = (0, A(i)) leads to the induced current J^uced = (^j ~§ A), 
which describes global color oscillations in the plasma. A more general Ansatz, A h 



A^(p ■ x — oot), has been presented in [13| to describe non-Abelian plane waves in the plasma. 
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Figure 1: J{x) and K(x) obtained by integrating numerically equations (|4.9| ), starting with 
regular boundary conditions at infinity. The plain and dashed lines represent different rates 
of approach to the (Yang-Mills vacuum) asymptotes at x = oo, with K = 1. Analysis of the 
oscillations in K near the origin [see the plots of K(x) / y/x vs. ln(x)] verifies the analytical 
form of the aymptotes (|4.10 ) . 



5. QGP Response and T = Topology 

In order to solve completely the time-depependent non-Abelian Kubo equation, one may 
need new insight, a possible source of which is provided by zero-temperature, topological gauge 
theory. Indeed, the gauge invariance condition for the generating functional of hard thermal 
loops in QCD 4 can be identified with the equation of motion for the gauge field in T = 0, d = 3 
topological Chern-Simons theory [H. 

Let us see how this works in some details. We first perform a gauge transformation 
S g A^ = d^oj + [A^jLu], with parameter to = —it a u a , of the HTLs generating functional (|3.1|) . 
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Imposing gauge invariance of ( |3.1|) , SgT^ThiA) = 0, implies 

J dfiq 5 g W(A + ) = 4tt y d A x (d A%)uj a . 



(5.1) 



Expressing 5 g W(A + ), upon integrating by parts, as 



S g W(A + ) 



d x 



5W(A A 
5A+ 



8 g A, 



dx <d- 



6W 
6A> 



+ 



.4 



+) 



SW 

Xa~ 



leads, upon introducing the combination / = £q- + A + , to the gauge invariance condition for 
hard thermal loops H, written now as: 



d f r, *i dA + 







(5.2) 



in the new space-time coordinates (u, v, x T ) defined by u = Qtx^, v = Q+x^, and Q + • x T = 0. 

To establish the relation with a topological action, one proceeds following ||. First, 
perform a Wick rotation into Euclidean 1R 4 space-time, and rename the coordinates u = z, 
v = z. Next, introduce the notation A + = a z and define a 2 = —f. Last, enforce the axial gauge 
condition a = 0. The gauge invariance condition for HTLs ( |5.2j ) can then be rewritten as 



d s a z - d z a- z + [a 2 , a z ] = . 



(5.3) 



This is the zero-curvature condition F zz = for the gauge field strength tensor F zz . I now 
switch to the Landau gauge, d^A^ = 0, for purposes of illustration. In that gauge, the zero- 
curvature condition ( |5.3| ) takes the form F^ u = d fl A u — d u A^ + [A^, A v ] = 0. This is just the 
equation of motion for the gauge field A^ with the classical action 



S 



dx e 



iivp 



A a r) A a A A a A b A c f° 



(5.4) 



The latter action is the one of topological Chern-Simons theory, a zero-temperature gauge field 
theory in three space-time dimensions. "Topological" refers to the fact that the Lorentz indices 
in ( |5.4| ) are saturated without using the metric tensor g^ u . As a consequence, the energy- 
momentum tensor, computed as = j^-, doesn't get any contributions from the action 
(|5.4|). Only the (non-topological, i.e., metric-dependent) gauge fixing and ghost terms that one 
adds to the action in order to quantize the theory contribute to T^, but such contributions do 
not represent physical energy-momentum. 

Chern-Simons theory, together with other topological field theories, has been intensively 
studied in recent years JI4"| , both due to its mathematical and physical interests. Chern-Simons 
theory describes anyons (particles with non-integer spins) and is a starting point for modeling 
the quantum Hall effect and high-T c superconductivity. Chern-Simons theory is known to 
possess the appealing property of perturbative finiteness [TH| , which in turn is connected to an 



interesting (twisted) form of supersymmetry [|16|, [15[ that is peculiar to topological theories. 
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The question wether such properties could be of interest in the realm of high-T QCD4 is an 
open, and interesting one. In particular, the twisted topological supersymmetry of Chern- 
Simons theory |17j could reveal a useful tool to shed new light on the relationship between high 



temperature in field theory and (untwisted) physical supersymmetry. 
6. Conclusions 

There are indications that the quark-gluon plasma may have been observed experimentally 
in heavy- ions collisions. However, a satisfactory theory of the QGP phase is still lacking. The 
best available theoretical framework is high-temperature QCD, in which hard thermal effects 
are dominant. 

The response of a quark-gluon plasma to a hard thermal perturbations is described by 
the non-Abelian Kubo equation. The static case is completely solved. In particular, there do 
not exist hard thermal solitons. This general result is confirmed by numerical computations, 
within a radial SU (2) Ansatz. The time-dependent case is less well understood. 

New insight into this field may hence be needed, and could indeed be provided by the 
close relationship that exists between hard thermal QCD 4 and three-dimensional Chern-Simons 
theory, a topological gauge theory at zero temperature. Chern-Simons theory is perturbatively 
finite. Furthermore, its gauge-fixed action exhibits a twisted supersymmetry which relates 
"bosons" (fields of even ghost number) to "fermions" (fields of odd ghost number), in very close 
analogy to usual (untwisted) supersymmetry. This leads us to the following, open questions. 
(1) Can the twisted supersymmetry of Chern-Simons theory be translated into the language of 
high-temperature QCD, and what does it mean in that context? (2) Can this shed new light 
onto the interplay between supersymmetry and high-temperature in field theory? 

Acknowledgements: The author expresses his gratitude to R. Jackiw and C. Manuel for 
interesting discussions and comments on the manuscript. 
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